Comment on "Foundation of Statistical Mechanics under Experimentally Realistic 

Conditions" 
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In a recent influential letter [lj , Reimann demonstrates 
equilibration of isolated quantum many-body systems 
from non-equilibrium initial states by proving the fol- 
lowing result: 


o\ ee [Tr{p(t)A} 
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where p m „(0) = (m\p(0)\n) denotes the initial den- 
sity matrix of the system under the eigenbasis {\n}} 
of the Hamiltonian H with eigenenergy {E n }, p(t) — 
E m „Pmn(0)e l ( E "- £ ™) t |m)(n| (setting h = 1) is the 
density matrix at time t, 777 denotes the time average, 
Peq = Pnn(0)\n)(n\ , and Aa is the working range of 
the observable A By assuming the conditions that 
1) the system's initial state has a broadly spread energy 
level population such that the quantity J2 n PnnO-0 * s ex ~ 
ponentially small with the number of degrees of freedom, 
and 2) the observable A has a bounded A^, the author 
shows that Eq. (1) implies that the time averaged vari- 
ance of the observable A from its equilibrium value is ex- 
ponentially small. This is a strong result. Unfortunately, 
we show here that the author's proof of the major result 
given by Eq. (1) is flawed and invalid for generic many- 
body systems under experimentally realistic conditions. 

There is a critical loophole in the proof of Eq. Q] 
in Ref. [1]. In simplifying the summation <j\ — 

Ek ra , n 4p* J ^n(0)e <[Bi - Bk+£ '"- B " li Eq.(ll) 
°f 111 (S means j k and m / n in the summation), 
the author drops all the terms except for the ones with 
j = n and k — m by exploiting a questionable "fact" 
that e l l E j- E k+E m -E n ]t van i snes f or a \\ the other terms 
under generic Hamiltonians. This "fact" is not valid for 
many-body systems under experimental realistic condi- 
tions. To see this clearly, let us first look at systems with 
quadratic Hamiltonians, which represent an important 
class of physical systems where the study of equilibration 
has recently attracted great theoretical and experimen- 
tal interest (see Refs. [2-4] in Reimann's paper). The 
quadratic systems can always be diagonalized as a set 
of non-interacting quasi-particle modes. The system has 
eigenenergy E n — J2i n i^i, where rij is the quasi-particle 
number of the i th mode with frequency Wj. The condi- 
tion Ej - Ek + E m - E„ = J2iti* - + m i - n i)ui = 
does not necessarily require ji = rii and fc, = m;. For in- 
stance, Ej — Ek + E m — E n — can be easily satisfied with 
ji+rrii — ki+rii (ji, rrii, fc,, n, are non-negative integers), 
which already leads to an exponentially large number 
of other solutions with non-vanishing e l ^ E ^ Eh+Em ~ En ^ t . 
Keeping only the terms with j = n and k — m in the 
summation drops out many (exponentially large in the 
number of modes) other terms, which leads to the artifact 


that a\ is bounded by an exponentially small quantity. 

One may argue that the special spectrum of the 
quadratic Hamiltonians leads to the above problem and 
hope that the terms with (j, k) ^ (n, m) can still be 
dropped for generic non-integrable Hamiltonians. This 
is not true and we show here that dropping the terms 
with (J, k) 7^ (n, m) is unjustified for any many-body 
Hamiltonians under experimentally realistic conditions. 
Note that for a large system with iV particles, typically 
its energy spread AE scales up linearly with the parti- 
cle number N, but the number of energy levels scales 
up exponentially with N. So, for many-body system, its 
energy spectrum is always exponentially dense, with the 
energy difference between adjacent levels going down ex- 
ponentially to zero. For any systems of physical rele- 
vance, the time window T for doing the average over 
e i[Ej-E k +E m -E n ]t j s finite. For instance, T cannot ex- 
ceed the age of the universe. However, because of the 
exponentially dense energy spectrum of many-body sys- 
tems, for any given small energy interval e ~ 1/T, there 
are always exponentially large number of near resonant 
terms with \Ej — Ek + E rn — E n \ < e which give non- 
vanishing e l [ E i- E k+E m -E n ]t ^ an( j a rj these terms cannot 
be dropped in the summation. In other words, to resolve 
the energy difference between adjacent levels for many- 
body systems, even for just hundreds of particles, the 
time window T required for the average 777 needs to be 
far beyond the age of the universe, which is clearly of no 
physical relevance. So the claim in Ref. [l[ for derivation 
of equilibrium under experimentally realistic conditions 
is unjustified. 

Although there are many recent works showing either 
numerically or experimentally that some isolated quan- 
tum systems appear to relax to p nn (0), much more com- 
plicated requirements on initial state and/or Hamiltonian 
are usually required for relaxation. The two conditions 
in Ref. [l[ are too simple and insufficient to guarantee re- 
laxation due to the aforementioned loophole in the proof. 

This work was supported by the NBRPC (973 Pro- 
gram) 2011CBA00300 (2011CBA00302), the DARPA 
OLE program, the IARPA MUSIQC program, the ARO 
and the AFOSR MURI program. 

Zhe-Xuan Gong and L. -M. Duan 
Department of Physics and MCTP, University of Michi- 
gan, Ann Arbor, Michigan 48109 and 
Center for Quantum Information, HIS, Tsinghua Univer- 
sity, Beijing, China 


[1] P. Reimann, Phys. Rev. Lett. 101, 190403 (2008). 


